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Excited electrons in the conduction band of germanium collect into four energy minima, or valleys, 
in momentum space. These local minima have highly anisotropic mass tensors which cause the elec¬ 
trons to travel in directions which are oblique to an applied electric field at sub-Kelvin temperatures 
and low electric fields, in contrast to the more isotropic behavior of the holes. This experiment pro¬ 
duces, for the first time, a full two-dimensional image of the oblique electron and hole propagation 
and the quantum transitions of electrons between valleys for electric fields oriented along the [0,0,1] 
direction. Charge carriers are excited with a focused laser pulse on one face of a germanium crystal 
and then drifted through the crystal by a uniform electric field of strength between 0.5 and 6 V/cm. 

The pattern of charge density arriving on the opposite face is used to reconstruct the trajectories of 
the carriers. Measurements of the two-dimensional pattern of charge density are compared in detail 
with Monte Carlo simulations developed for the Cryogenic Dark Matter Search (SuperCDMS) to 
model the transport of charge carriers in high-purity germanium detectors. 

PACS numbers: 72.10.Di, 72.20.Jv, 72.20.-i, 72.20.Dp, 72.80.Cw, 71.18.+y 


INTRODUCTION 


Because germanium is an indirect-gap semiconductor, 
the energy minimum of its conduction band does not oc¬ 
cur at zero momentum. Rather, the conduction band 
has four minima, or valleys, located at the edges of the 
Brillouin zone in the four (1,1,1) crystal directions. The 
effective masses of the electrons in these valleys are highly 
anisotropic, and take on values of 1.58 and 0.081 times 
the mass of the electron in vacuum in the directions par¬ 
allel and perpendicular to the (1,1,1) directions, respec¬ 
tively. Ilia Because of the anisotropic mass, each elec¬ 
tron’s acceleration is rarely parallel to the applied force; 
the direction of acceleration depends on the mass tensor 
of the valley the electron is occupying. As a result, an 
initially localized group of electrons in a uniform elec¬ 
tric field will spatially separate into four clusters, each 
consisting of electrons occupying one of the four valleys. 

At high temperatures, or in the presence of high elec¬ 
tric fields or impurity concentrations, electrons will un¬ 
dergo frequent quantum transitions between these four 
valleys, resulting in a nearly isotropic electron mobility. 
However, if the temperature, electric field, and impurity 
concentrations are all sufficiently low, the mean-free-path 
for inter-valley scattering can exceed the sample size, and 
the mass anisotropy will become apparent. This effect 
was the cause of a surprising asymmetry observed be¬ 
tween the collection of electrons and holes in the ger¬ 
manium crystals used as dark matter detectors in the 
Cryogenic Dark Matter Search (CDMS). [3] To model 
this asymmetry, we developed Monte Carlo simulations of 


charge transport in our germanium crystals. The follow¬ 
ing experiment was designed to verify these simulations 
in detail by directly imaging the propagation of charge 
carriers. 



(a) Illustration of scanning system. 



(b) Photograph of a raster scan. 


FIG. 1. Laser Scanning System: (a) Schematic of the 
crystal sample and laser scanning system, (b) Benchtop pho¬ 
tograph of the scanning system with a raster-scan across the 
surface of the crystal and the copper-clad GIO crystal holder, 
shown here outside of the He-3 cryostat. Gopper shielding 
removed for photograph. 

EXPERIMENTAL SETUP 

The crystal under test was cut from a 3.89mm thick 
wafer of high purity germanium. [4] The front and back 
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faces are Icmxlcm, and lie in the (0,0,1) crystal plane, 
while the sides lie in the (1,1,0) and (1,1,0) planes. The 
illuminated (front) face of the crystal is patterned with 
an aluminum-tungsten mesh electrode, with 10/im pitch 
and 20% coverage. [5] The non-illuminated (back) face is 
covered by two solid electrodes, one of which is circular 
with a diameter of 160/im, and separated by a 10/im gap 
from the other electrode, which covers the rest of the face 
of the crystal. Both electrodes on the non-illuminated 
face are connected to ground through 10.8MO resistors. 

A benchtop photograph of the crystal sample is shown 
in Fig. [^. The sample is mounted to a copper baseplate 
(bottom of photograph), and is cooled to 600mK under 
vacuum in a He-3 cryostat. The crystal sample is shielded 
from both electrical noise and 4K thermal radiation by a 
copper enclosure mounted on the baseplate. 

To generate free charge carriers, the illuminated face 
is exposed to a 20/iW, 100ns pulse of 650nm laser light, 
focused to a 60/im diameter spot, which creates a cloud 
of electron-hole pairs near the illuminated surface. Car¬ 
riers are propagated through the bulk of the crystal in 
the uniform electric field induced by the DC bias voltage 
of the mesh electrode on the illuminated face. The sign 
of the bias voltage determines whether electrons or holes 
propagate through the bulk of the crystal. Once col¬ 
lected by the electrodes on the non-illuminated face, the 
charge carriers produce proportional pulses in voltage, 
which are measured by high-impedance MESFET ampli¬ 
fiers, mounted on the 4K stage in the He-3 cryostat. EE] 

The two-dimensional pattern of charge density is de¬ 
termined by measuring the charge collected by the small, 
circular electrode as a function of the position of the exci¬ 
tation point on the illuminated face. Because translation 
of the excitation point causes a corresponding translation 
of the charge density pattern, this procedure is equiva¬ 
lent to keeping the excitation point fixed, while moving 
the small, circular electrode to probe different regions 
of the charge density pattern. The position of the exci¬ 
tation point is controlled by means of a Micro Electro- 
Mechanical System (MEMS) mirror from Mirrorcle Tech¬ 
nologies, Inc. [8] The MEMS mirror can be tilted along 
two axes by computer control, and was modified by the 
vendor for operation below IK. 

RESULTS 

The experimental results are shown in Eigs. and 
[^, which plot the voltage amplitude of the small circu¬ 
lar electrode as a function of the position of the point 
of laser excitation. These data provide a map of the 
two-dimensional pattern of charge arriving on the non- 
illuminated face of the crystal, and are compared to the 
simulated distribution of charge carriers calculated by 
Monte Carlo simulation, shown in Eigs. It and Ht- 

The comparison between the Monte Carlo simulation 
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(c) Data (solid) vs. Simulation (dotted) 

FIG. 2. Hole Charge Density Patterns: (a): Data, (b): 
Simulation, (c): One-dimensional projection of charge den¬ 
sity onto x-axis. The data (solid, blue) are compared to the 
simulation incorporating electrostatic repulsion (dotted, red), 
and the simulation with no repulsion (dotted, green). The 
horizontal scale ranges from -4mm to -|-4mm. The vertical 
scale is arbitrary. 


and the experimental data can be more clearly seen in 
Figs. an d|l, which show one-dimensional projections 
of the two-dimensional charge density patterns, for both 
the data (solid line) and the simulation (dotted line), in 
the presence of several different electric field strengths. 

The effects of the electron mass anisotropy are clearly 
evident in Fig. |§. At low electric fields (less than 
^3V/cm), the electrons separate into four groups, each 
corresponding to one of the conduction band minima, or 
valleys. At high electric fields (greater than ^4V/cm), 
the inter-valley scattering induced by phonon emission 
washes out the effects of the anisotropic mass, and the 
charge density pattern approaches a Gaussian distribu¬ 
tion. 

Figure shows the mean path length between electron 
inter-valley scattering events, projected along the [0, 0,1] 
direction, as a function of electric field strength. The 
data from our experiment are shown in blue, and are 
in good agreement with the Redl simulation, shown in 
green. Note that the parameters which determine the 
electron scattering rate in the Redl simulation were taken 
from reference m, and were not modified to fit our data. 
For comparison, we also show data from references m 
and m, plotted on the same axes. 
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(a) Electron Data 




(c) Data (solid blue) vs. Simulation (dotted red) 


FIG. 3. Electron Charge Density Patterns: (a): Data, 
(b): Redl simulation, (c): One-dimensional projection of 
charge density onto a diagonal axis. The data (solid, blue) are 
compared to the Redl simulation employing the Herring-Vogt 
approximation (dotted, red). The horizontal scale ranges 
from -4mm to +4mm. The vertical scale is arbitrary. 


MONTE CARLO SIMULATION 


To obtain good agreement between the Monte Carlo 
simulation by Peter Redl and the experimental data, sev¬ 
eral effects need to be taken into account. The first ef¬ 
fect is a distortion in the locations of the centers of the 
four electron clusters. This distortion is explained by a 
2 degree offset between the crystal axes and the faces 
of the crystal, which is consistent with the tolerances of 
our crystal fabrication process and with subsequent X- 
ray diffraction measurements of the crystal axes. 

The second effect is the lateral spreading of the carriers 
due to electrostatic repulsion. Without accounting for 
this effect, our original Monte Carlo simulation for the 
holes produced charge distributions that were too nar¬ 
row (see Fig. [^). The additional spreading seen in our 
data can be reproduced in our simulations by incorpo¬ 
rating the effect of electrostatic repulsion. Based on the 
total charge collected by the large electrode on the non- 
illuminated face, we estimate 1.9±.4x 10^ charge carriers 
are collected per laser pulse. [9] To obtain good agreement 
with the data, the simulation contains 2 ± 1 x 10^ holes, 
assuming the initial charge cloud diameter after the sep¬ 
aration of opposite charges is equal in size to the laser 
spot (60/im). 

However, even after accounting for electrostatic repul- 


Inter-Valley Scattering Length vs. Electric Field 



FIG. 4. Electron inter-valley scattering length, pro¬ 
jected along [0,0,1], as a function of electric field strength. 
The data (blue) are compared to the Redl simulation (green) 
and the Sundqvist simulation (black). Also included is the IV 
scattering length calculated using data from references m 
and [T^. 



(a) Redl Simulation (b) Sundqvist (c) Sundqvist Sim. 

Simulation with Repulsion 


FIG. 5. Comparison of Redl and Sundqvist Simula¬ 
tions: (a), (b), and (c) show scatter plots of the hnal x and y 
electron locations in three simulations. The axes extend from 
-4mm to +4mm in the x and y-directions. (a) Reid sim¬ 
ulation using Herring-Vogt approximation, which does not 
conserve momentum, (b) Sundqvist simulation, which cor¬ 
rectly conserves momentum during phonon scattering. [HHH] 
(c) Sundqvist simulation with correction for electrostatic re¬ 
pulsion. For all simulations, the electric field strength was 
2V/cm, and no impurity scattering was included. 


sion, a discrepancy remains in the case of the electrons, 
as the electron charge distribution is narrower in the data 
than in the Redl simulation (compare Figs. &[^). We 
attribute this discrepancy to the fact that the Redl simu¬ 
lation uses the Herring-Vogt approximation. [TOl [13] This 
approximation provides a 20x faster means of simulating 
the phonon emission process at the expense of violating 
momentum conservation. This non-conservation of mo¬ 
mentum is likely the cause of the extra lateral diffusion 
seen in the simulated charge density patterns. 

To estimate the systematic error introduced by the 
Her ring-Vogt approximation, we use a simulation writ- 
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FIG. 6. Comparison of Lateral Spreading of Electron 
Charge Density: (a) and (b) show plots of one-dimensional 
charge density of a single electron cluster, projected along the 
narrow and wide principle axes, respectively. A Gaussian fit 
to the data is shown in black, and is compared to three differ¬ 
ent simulations. For both simulations and data, the electric 
field strength was 2V/cm. 

ten by Kyle Sundavist pAUTB] which explicitly conserves 
crystal momentum. Figure shows a comparison of the 
Redl simulation (a), with the more accurate, but slower, 
Sundqvist simulation (b). The widths of the distribu¬ 
tions of electrons in a single cluster are shown in Fig. 
for both simulations, and for the data. 

Using the Sundqvist simulation as a standard of com¬ 
parison, we are able to observe and model the extra lat¬ 
eral spreading of the electrons due to electrostatic repul¬ 
sion. The repulsion causes the centers of the four electron 
clusters to spread further from each other, and causes the 
clusters to spread more in the low-mass directions. Fig¬ 
ure If shows an approximate correction to the Sundqvist 
simulation for the effects of electrostatic repulsion, gener¬ 
ated by post-processing the electron trajectories. Figure 

shows how this effect brings the Sundqvist simulation 
into agreement with the experimental data, assuming the 
same number of charge-carriers per pulse (2 ± 1 x 10^) as 
we assume for holes. 

We also note that, unlike the Redl simulation, the 
Sundqvist simulation does not include the effects of im¬ 
purity scattering. This fact is apparent in Fig. where 
the IV scattering length for the Sundqvist simulation di¬ 
verges for low electric field. However, a discrepancy still 
exists between the Redl and Sundqvist simulations even 
at high fields where scattering is dominated by phonon 
emission. This discrepancy remains unexplained. 


CONCLUSIONS 

Using our novel cryogenic scanning apparatus, we 
have, for the first time, imaged the two-dimensional 
charge density distributions of both electrons and holes 
in high purity germanium at sub-Kelvin temperatures. 
These results provide a useful check on the accuracy of 
our Monte Carlo simulations and provide motivation for 
improvements. Plans for future work include investiga¬ 
tion of the scattering and trapping rates of charge car¬ 
riers in multiple germanium and silicon crystals, which 
will improve our understanding of the low-temperature 
condensed-matter physics relevant to the design and 
analysis of the germanium and silicon detectors used by 
the CDMS collaboration for dark matter searches. 
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